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(54) System and method for predicting the timing of future service events of a product 



(57) A system (28) and method for predicting timing 
of future service events of a product. A database (30) 
contains a plurality of service information and perform- 
ance information for the product. A statistical analyzer 
(34) analyzes the plurality of processed service informa- 
tion to determine a plurality of compartment failure in- 



formation. A performance deterioration rate analyzer 
(42) analyzes the performance deterioration rate of the 
product from the plurality of service information and per- 
formance information. A simulator (44), simulates a dis- 
tribution of future service events of the.product accord- 
ing to the plurality of compartment failure information 
and the performance deterioration rate analysis. 
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Description 

[0001] This disclosure relates generally to servicing products and systems and more particularly to predicting the 
timing of future service events of a product or a system. 

5 [0002] The market for long-term contractual agreements has grown at high rates over recent years for many of 
today's service organizations. As the service organizations establish long-term contractual agreements with their cus- 
tomers, it becomes important to understand the expected costs and risks associated with the pricing of service contracts 
and portfolio management of the contracts. In addition, the service organizations need to have an understanding of 
the planning of repairs (shop workload planning) and how the introduction of new technology will affect their service 

10 contracts. In order to analyze these issues, it is necessary to correctly model the underlying behavior of the product 
or system so that each can be serviced in the most cost-effective manner. 

[0003] Cu rrently available analytical practices are unable to accurately model service requirements for complex prod- 
ucts or systems. Typically, these models contain poor cost information which result in the service organization ineffi- 
ciently managing the risk associated with their service portfolios, failing to respond to customer needs and new tech- 

15 nology, which all lead to lower long-term contract profitability. A standard time-series method is one particular approach 
that has been used to model the service requirements of repairable systems such as aircraft engines, automobiles, 
locomotives and other high tech products. This time-series method examines historical data obtained over a five to 
ten year period and forms a trend line on either system costs and/or number of repairs made to the system. The trend 
line is then used to predict future costs and number of repairs. A limitation with this time series method is that it does 

20 not give details of failures at a compartmental level. A compartment is a physical or performance related sub-system 
of the repairable product, which when it fails causes the product to require maintenance or servicing. Other limitations 
with the standard time series method is that it does not account for the life cycle of the repairable product and thus 
does not provide a distribution of the expected service events for the product. An analysis based on engineering rela- 
tionships to determine compartment parameters is another method used to model the service requirements of repair- 

25 able systems. A limitation with this analysis is that it is not well based in underlying statistics, and thus cannot be shown 
to accurately model the repairable product on an ongoing basis. 

[0004] In order to overcome the above problems, there is a need for an approach that can model the service require- 
ments of repairable systems that is accurate and has a comprehensive statistical framework. Such an approach will 
lead to better cost projections, more realistic and effective risk management, new technology introduction and day-to- 

30 day service that is more responsive to customer needs and higher long-term contract profitability. 

[0005] In accordance with one embodiment of this disclosure, there is a system for predicting the timing of a future 
service event of a product formed from a plurality of compartments. The system comprises a database that contains 
a plurality of service information and a plurality of performance information for the product. A statistical analyzer ana- 
. lyzes the plurality of service information to determine a plurality of compartment failure information. A performance 

35 deterioration rate analyzer analyzes the performance deterioration rate of the product from the plurality of service 
information and performance information. A simulator, simulates a distribution of future service events of the product 
according to the plurality of compartment failure information and the performance deterioration rate analysis. 
[0006] Similarly, in this disclosure there is a method for predicting the timing of a future service event of a product 
formed from a plurality of compartments. The method comprises storing a plurality of service information and a plurality 

40 of performance information for the product; analyzing the plurality of service information to determine a plurality of 
compartment failure information; performing a deterioration rate analysis of the product from the plurality of service 
information and performance information; and simulating a distribution of future service events of the product according 
to the plurality of compartment failure information and the deterioration rate analysis. 

[0007] Also, in this disclosure there is a computer-readable medium storing computer instructions for instructing a 
45 computer to predict the timing of a future service event of a product formed from a plurality of compartments. The 
computer instructions comprise storing a plurality of service information and a plurality of performance information for 
the product; analyzing the plurality of service information to determine a plurality of compartment failure information; 
performing a deterioration rate analysis of the product from the plurality of service information and performance infor- 
mation and simulating a distribution of future service events of the product according to the plurality of compartment 
50 failure information and the deterioration rate analysis. 

[0008] An embodiment of the invention will now be described, by way of example, with reference to the accompanying 
drawings, in which: 

Fig. 1 shows a schematic of a general purpose computer system in which a system for predicting the timing of 
55 future service events of a product operates; 

Fig. 2 shows a schematic diagram of a system for predicting the timing of future service events of a product that 
operates on the computer system shown in Fig. 1 ; 
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Fig. 3 shows a flow chart describing actions performed by the system shown in Fig. 2; 

Fig. 4 shows a flow chart describing the actions performed by the statistical analyzer shown in Fig. 2; 

5 Fig. 5 shows a flow chart describing the actions performed by the performance deterioration rate analyzer shown 

in Fig. 2; 

Figs. 6a-6b show examples of plots that describe some of the performance information stored in the performance 
historical database shown in Fig. 2; 

10 

Fig. 7 shows a flow chart describing the actions performed by the simulator shown in Fig. 2; and 

Fig. 8 shows a flow diagram describing the validating actions performed by the system shown in Fig. 2. 

15 [0009] Fig. 1 shows a schematic of a general-purpose computer system 1 0 in which a system for predicting the timing 
of future service events of a product operates. The computer system 1 0 generally comprises a processor 1 2, a memory 
14, input/output devices, and data pathways (e.g., buses) 16 connecting the processor, memory and input/output de- 
vices. The processor 12 accepts instructions and data from the memory 14 and performs various calculations. The 
processor 12 includes an arithmetic logic unit (ALU) that performs arithmetic and logical operations and a control unit 

20 that extracts instructions from memory 14 and decodes and executes them, calling on the ALU when necessary. The 
memory 14 generally includes a random-access memory (RAM) and a read-only memory (ROM), however, there may 
be other types of memory such as programmable read-only memory (PROM), erasable programmable read-only mem- 
ory (EPROM) and electrically erasable programmable read-only memory (EEPROM). Also, the memory 14 preferably 
contains an operating system, which executes on the processor 12. The operating system performs basic tasks that 

25 include recognizing input, sending output to output devices, keeping track of files and directories and controlling various 
peripheral devices. 

[0010] The input/output devices comprise a keyboard 18 and a mouse 20 that enter data and instructions into the 
computer system 10. A display 22 allows a user to see what the computer has accomplished. Other output devices 
could include a printer, plotter, synthesizer and speakers. A modem or network card 24 enables the computer system 

30 1 o to access other computers and resources on a network. A mass storage device 26 allows the computer system 10 
to permanently retain large amounts of data. The mass storage device may include all types of disk drives such as 
floppy disks, hard disks and optical disks, as well as tape drives that can read and write data onto a tape that could 
include digital audio tapes (DAT), digital linear tapes (DLT), or other magnetically coded media. The above-described 
computer system 1 0 can take the form of a hand-held digital computer, personal digital assistant computer, personal 

35 computer, workstation, mini-computer, mainframe computer and supercomputer. 

[0011] Fig. 2 shows a schematic diagram of a system 28 for predicting the timing of future service events of a product 
that operates on the computer system 1 0 shown in Fig. 1 . In the system 28, a service database 30 stores a plurality 
of service information for the product. The plurality of service information varies depending on the product. Generally, 
the plurality of service information will comprise information such as compartment definitions of the product (i.e., a 

40 physical or performance related subsystem that is considered as a unit, which when it fails requires that the product 
needs maintenance or servicing), repair history of the product (e.g., dates of service events, types of service events, 
time for a compartment to fail, etc.), as well as any factors which may play a role in explaining the length of time which 
passes between service events (e.g., environment, operating conditions of the product, product configuration, equip- 
ment age, etc.). Other factors may include types of maintenance, cycle time of the product, usage of the product, 

45 contract terms and conditions, equipment age and vintage, etc. In this disclosure, the product will be described with 
reference to an aircraft engine; however, other products such as a power system, a locomotive or any other electrical, 
chemical or mechanical products, where it is desirable to predict the timing of future service events, may be used. 
[0012] Referring to Fig. 2, a preprocessor 32 processes the plurality of service information into a predetermined 
format. The preprocessing includes extracting the plurality of service information from the service database 30; as- 

50 signing each data record in the database to a compartment depending upon the engineering labeled removal cause; 
creating new variables from existing variables (e.g., censoring variables, customer indicator variables), deleting outliers 
and producing summary statistics of the data set (e.g., number of records for each compartment). After performing 
these acts, the preprocessor 32 generates a plurality of data files according to the service information, wherein each 
of these data files are formatted as SAS data sets. 

55 [0013] A statistical analyzer 34 analyzes the plurality of processed service information to determine a plurality of 
compartment failure information. The compartment failure information may include statistically significant compartment 
failure variables and the associated compartment time-to-failure coefficients. Compartment failure variables are vari- 
ables that affect the time between service or maintenance events. For example, in the aircraft engine scenario, the 
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thrust rating of the engine and the environment that the engine flies in are examples of possible statistically significant 
compartment failure variables. Compartment time-to-failure coefficients are coefficients that are applied to each of the 
compartment failure variables. The compartment failure variables and associated compartment time-to-failure coeffi- 
cients are used to determine the time between servicing events for the compartments. In addition, the statistical analyzer 
5 34 uses this information to determine which compartment failure variables influence service events and estimate time- 
to-failure distributions for the compartments. 

[0014] The statistical analyzer 34 comprises several scripts that enable it to perform the aforementioned functions 
as well as some additional functions. One particular script that the statistical analyzer 34 uses is a service analysis 
script that executes a plurality of statistical procedures. The plurality of statistical procedures may comprise a multi- 

10 variate regression and/or a correlation analysis. Both the multi-variate regression and correlation can determine which 
compartment failure variables influence the timing of service events and estimate the time-to-failure distributions for 
the compartments. The statistical analyzer 34 uses other scripts to output this information as well as statistical diag- 
nostics 36 and residual plots 38. The output from the multi-variate regression and/or correlation analysis may include 
the compartment time-to-failure coefficients for each compartment associated with the product. Other standard output 

15 from the multi-variate regression and/or correlation analysis may include a standard error associated with each of the 
compartment time-to-failure coefficients and a P value, which is an indication of whether a particular variable has a 
significant effect on the time between occurrences of service events. 

[0015] Statistical diagnostics that may be outputted include goodness-of-fit metrics and collinearity diagnostics. 
These enable a user to determine the most appropriate compartment model in the statistical analyzer 34. Residual 
20 plots enable a user of the system 28 to determine how well the regression model fits the service information data. 
Generally, the residual plots are defined as the difference between the actual time until servicing values and the pre- 
dicted time until servicing values. A small residual value is an indication that the regression or correlation analysis was 
a good fit, whereas a large residual value is an indication that the fit may be improved. 

[0016] In addition to residual plots, the statistical analyzer 34 may use another script to output information such as 
25 probability plots, which enables one to assess whether the assumed life distribution for each compartment is appropriate 
or not. Another script may be used to generate a plot of residuals versus each variable that affects the time that a 
service event occurs. 

[0017] Referring again to Fig. 2, the system 28 also comprises a performance historical database 40 that includes 
a plurality of performance information obtained from the product while in operation. As mentioned above, this disclosure 

so is described with reference to an aircraft engine. Therefore, the performance information can be acquired by using any 
of a plurality of data acquisition devices such as sensors and transducers. After the data have been obtained, the data 
acquisition devices can transfer the data to a remote monitoring facility for storage and evaluation. Also, it is possible 
to have the data from the data acquisition devices manually recorded and entered into the performance historical 
database 40. The plurality of performance information includes but is not limited to information such as performance 

35 characteristic values (e.g., exhaust gas temperature, EGT), initial data levels after servicing, current data levels, dates 
at which the product is serviced, and variables that affect the servicing of a subset of compartments of the product. 
For this aircraft engine embodiment, the variables may include flight leg, engine thrust, customer, engine model, engine 
series, etc. All of this performance information is described below in more detail. 

[0018] A performance deterioration rate analyzer 42 analyzes the performance characteristic values of the product 

40 from both the plurality of service information and performance information. The performance deterioration rate analyzer 
42 comprises a statistical analysis script that relates a subset of compartments of the product according to time. For 
purposes of describing the performance deterioration rate analyzer 42, time is the amount of time that the aircraft 
engine is in use. Time can be measured by variables such as cycles or hours. The statistical analysis script generates 
an estimated deterioration rate curve for a subset of compartments of the product. The performance deterioration rate 

^5 analyzer 42 further comprises a transformer that transforms each estimated deterioration rate curve for a compartment 
to a performance life distribution. A performance life distribution is a statistical distribution representing the statistical 
properties of the time between servicing events and is estimated using performance data as opposed to service data. 
The performance life distribution is in the same form as the estimated time-to-failure distributions for the compartments 
determined by the statistical analyzer 34. 

so [0019] A simulator 44 simulates the future service events of the product according to the plurality of compartment 
failure information generated by the statistical analyzer 34 and the performance life distribution generated by the per- 
formance deterioration rate analyzer 42. The simulation results in a forecast of the timing of the future service events. 
In particular, the simulator 44 takes the compartment time-to-failure coefficients from the statistical analyzer 34 and 
determines a Weibull distribution for each compartment defined for the product. Also, the simulator 44 takes the per- 

55 formance life distribution from the performance deterioration rate analyzer 42 and determines a Weibull distribution for 
each of the associated compartments defined for the product. The simulator 44 then uses the compartment distributions 
to determine the overall distribution for the product. Generally, the simulator uses a discrete event-driven Monte Carlo 
simulation to perform the above operations. After performing the simulation, the simulator 44 generates several outputs. 
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For instance, one output is the contract output 46, which typically comprises the following: maintenance event distri- 
bution parameters overtime; demand distributions spare/leased equipment; and equipment performance distributions 
(e.g., aircraft engine time on wing). The simulator 44 is not limited to these outputs and it is possible to use other 
outputs if desired. 

5 [0020] In order to understand the performance of the simulator 44, the system 28 uses a validator, which can be part 
of the statistical analyzer 34 or the simulator or separate from both. The validator contains a validation script prepared 
for a case study done for the product. For purposes of this disclosure, a case study is defined as any subset of historical 
service event data that is used for model validation. For example, the service events that took place on a group of 
randomly chosen systems over the past year may serve as a case study, and during validation, a comparison is made 

10 between the number of service events projected by the model for these systems over this period of time and the actual 
number of service events that took place. The validation script will compare the compartment distributions determined 
by the simulator 44 to the distributions that actually happened in the case study. After making the comparison, the 
validator generates a series of graphical outputs 48 on availability and reliability. In a preferred embodiment, three sets 
of reliability graphs are generated. The first set of reliability graphs are relative frequency histograms of the actual 

15 compartment distributions for each of the first four shop visits for the case study. Overlaid on these relative frequency 
histograms are the compartment distributions determined by the simulator 44. The second set of reliability graphs are 
relative frequency histograms of the actual system level distributions for each of the first four shop visits for the case 
study. Overlaid on each of these relative frequency histograms are the system distributions determined by the simulator 
44. The third set of reliability graphs are non-parametric Kaplan-Meier estimated survival curves determined from both 

20 the actual system level distribution and the system level distribution determined by the simulator 44 for each of the first 
four shop visits. From these outputs, a user can generate a service plan forecast for the product that comprises time 
for scheduling service events. 

[0021] Fig. 3 shows a flow chart describing actions performed by the system 28 shown in Fig. 2. At block 50, a 
plurality of service information and performance information for the product stored in the service database and the 

25 performance historical database, respectively, are obtained. The preprocessor preprocesses the plurality of service 
information into a predetermined format at 52. The statistical analyzer analyzes the plurality of processed service 
information to determine a plurality of compartment failure information at 54. In particular, the statistical analyzer de- 
termines both the compartment time-to-failure coefficients and the compartment failure variables using the aforemen- 
tioned statistical procedures. The statistical analyzer outputs the compartment time-to-failure coefficients and the com- 

30 partment failure variables to the simulator and generates various residual plots. 

[0022] At the same time the service information are being preprocessed and analyzed by the preprocessor and the 
statistical analyzer, the performance information are simultaneously evaluated by the performance deterioration rate 
analyzer. If desired, it is also possible to have the performance information preprocessed in a manner similar to the 
service information. Regardless of whether the performance information are preprocessed, the performance deterio- 

35 ration rate analyzer runs a deterioration rate analysis at 56. As mentioned above, the deterioration rate analysis gen- 
erates an estimated deterioration rate curve for a subset of compartments of the product and transforms each estimated 
deterioration rate curve to a performance life distribution. 

[0023] After analyzing the service information and the performance information, the simulator simulates the future 
service events of the product according to the compartment failure information and the performance life distribution at 

40 58. In addition, the simulator forecasts or predicts the timing of the future service events at 60. As mentioned above, 
this information is in the form of distributions for each compartment that makes up the product. The validator compares 
the compartment distributions determined by the simulatorto the distributions that actually happened in the case study 
at 62. After making the comparison, the validator generates a series of graphical outputs on availability and reliability. 
[0024] Fig. 4 shows a flow chart describing the actions performed by the statistical analyzer shown in Fig. 2. At block 

45 64, the statistical analyzer obtains the service information from the preprocessor. The statistical analyzerthen generates 
compartment definitions forthe service information at 66. At 68, the statistical analyzer determines compartment failure 
information such as the statistically significant compartment failure variables and their associated compartment time- 
to-failure coefficients using the aforementioned statistical procedures. The statistical analyzerthen applies the com- 
partment time-to-failure coefficients to the compartment failure variables at 70. At block 72, the statistical analyzer 

50 generates the various statistical diagnostics for each compartment associated with the product. At block 74, the sta- 
tistical analyzer generates residual plots and probability plots and other types of plots if desired. As mentioned above, 
the statistical analyzer can generate other information such as standard errors associated with each of the compartment 
time-to-failure coefficients and a P value. . 

[0025] Fig. 5 shows a flow chart describing the actions performed by the performance deterioration rate analyzer 42 
55 shown in Fig. 2. Before running the deterioration rate analysis, the performance deterioration rate analyzer first obtains 
the plurality of service information and performance information from the service database and the performance his- 
torical database, respectively, at 76. As mentioned earlier, the plurality of performance information includes information 
such as performance characteristic values, initial data levels after servicing, current data levels, dates at which the 
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product is serviced, and variables that affect the servicing of each compartment of the product. 
[0026] Figs. 6a-6b show examples describing some of the above-mentioned performance information. In particular, 
Fig. 6a shows an example of a deterioration rate curve for a compartment of an aircraft engine. In this example, the 
compartment is EGT, however, other compartments could be used. For instance, an illustrative but non-exhaustive list 

5 could include delta exhaust gas temperature (dEGT), which is the deviation from the baseline EGT, fuel flow (WF), 
core speed (N2), and EGT divergence, divEGT, which is the difference of the EGT between the raw EGT of the engine 
in question and the mean of raw EGT of all engines. The deterioration curve shows the degradation of the compartment 
over time. Eventually, after a period of time, the compartment reaches a level that the degradation is severe enough 
to warrant servicing. The initial data level performance parameter is the initial level of the compartment after being 

10 serviced. In Fig. 6a, the initial level is about 50°F, and over time the level of the EGT margin will degrade. The historical 
trending levels are shown in Fig. 6a as data points. Fig. 6b shows an example of the raw EGT levels. In particular, Fig. 
6b shows a plot illustrating the removal level or redline for the EGT compartment. The removal level indicates an 
absolute time that the compartment reaches a predetermined level that necessitates the removal of the aircraft engine 
for servicing. 

15 [0027] Referring back to Fig. 5, after the service information and the performance information are obtained, the 
performance deterioration rate analyzer executes the statistical analysis script that relates each compartment of the 
product according to time at 78. Preferably, the statistical analysis runs a multi-variate regression analysis for each 
compartment of the product to identify variables that influence the time between servicing events. An illustrative ex- 
ample of a multi-variate regression analysis using a Weibull distribution is presented. The time (specified either in 

20 engine flight hours or engine cycles) between service events is represented as Y. As mentioned earlier, thrust (X1 ) and 
flight leg (X2) are two compartment failure variables that might influence the time between service events, Y, The multi- 
variate Weibull regression model takes the form: 

ln(V) = a + p 1 X1 +p 2 X2 + oe, 

25 

where In is the natural log function, e is an error term which follows the smallest extreme value distribution, and a, 
p 2 , and a are compartment failure parameters to be estimated from the service data. 

[0028] For fixed values of thrust (X1 ) and flight leg (X2), a Weibull distribution representing time between servicing 
30 events may be determined. For example, if <x=8.9, p^-0. 00003, ^=0.75, and o=0.5 where X1=23500 and X2=1.8, 
the Weibull distribution of the time between servicing events for this compartment would have a location (or scale) 
parameter (i.e. , the 63.2 nd failure percentile) of exp(8.9-0.00003*23500+0.75*1 .8)= 1 3975 and a shape parameter eq ual 
to 1/g=1/0.50=2.0. Although a Weibull regression analysis is described, other statistical analyses such as multiple non- 
linear and loglinear analyses could also be used. 
35 [0029] The performance deterioration rate analyzer generates estimated deterioration rate curves for a subset of 
compartments of the product at 80. The estimated deterioration rate curves are determined using a multi-variate re- 
gression and/or correlation statistical analysis. Denote the performance characteristic (e.g., EGT margin) as Y 
[0030] One example of a multi-variate regression and/or correlation analysis is presented using Y and time (as meas- 
ured in cycles) in the following model: 
40 ■ 

Y = a + fi^Cycies + e 

where e follows a normal distribution and a and p 1 are parameters to be estimated from the performance and service 
45 information data. The estimate of p 1 in this example is an estimated rate of deterioration for the performance charac- 
teristic, Y. For example, if P^ 0.003 then the performance characteristic, Y, deteriorates at approximately 3 degrees 
per 1000 cycles. . 

[0031] Next, the performance deterioration rate analyzer transforms each estimated deterioration rate curve for the 
respective compartment to a performance life distribution at 82. The performance life distribution is characterized by 
so a location (or scale) parameter and a shape parameter. In order to perform the transformation, estimates of a, fa , and 
a performance characteristic limit value (i.e., a value at which the compartment requires servicing) which we denote 
by EGTL, are required. The location (or scale) parameter of the performance life distribution is obtained using the 
following formula: 

55 
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[0032] The shape parameter of the performance life distribution is obtained using the following formula: 



Shape = 



ln(-ln(l -0.825)) 



15 




In 



20 [0033] As an example, suppose the estimate of a is 0, the estimate of ^=0.003 and the performance characteristic 
value is 60. Using the formula above, the location parameter value is estimated to be 23498. The shape parameter is 
estimated to be 2.27. The performance life distributions for all compartments are then transferred to the simulator at 84. 
[0034] Fig. 7 shows a flow chart describing the actions performed by the simulator 44 shown in Fig. 2. As mentioned 
earlier, the simulator 44 is interested in determining the distribution of failures at the product's system level so that the 

25 timeliness of future service events can be predicted. The simulator 44 is able to determine the distribution of failures 
at the product's system level because of the information provided by the statistical analyzer 34 and the performance 
deterioration rate analyzer 42. The information (i.e., time-to-failure coefficients and compartment variables) provided 
by the statistical analyzer 34 facilitates an understanding of each of the compartments that make up the product's 
system level and their relationship with each other, while the performance life distributions provided by the performance 

30 deterioration rate analyzer gives more information about probable service requirements. The simulator 44 uses thjs 
information to examine the system or aggregate level and predict the overall timing of service events for the product. 
Referring back to Fig. 7, the actions performed by the simulator begin at block 86, where random failure times (i.e., 
service events) for each compartment distribution are generated. From the randomly generated failure times, the min- 
imum of these values is found at 88. The simulator designates t f as the minimum time, where / is the compartment 

35 associated with this time value. The simulator then records the minimum time f, as the next failure time (i.e., time for 
a service event) for the system level at 90. At 92, the simulator determines whether there are any more system level 
failures needed. If so, then blocks 86-90 are repeated a large number of times. Once all of the iterations have been 
performed, the simulator forms a system level distribution from the failure times at 94. At 96, the simulator generates 
the output tables and the input report, while a graphical output is generated at 98. 

40 [0035] Fig. 8 shows a flow diagram describing the validating actions performed by the system shown in Fig. 1. In 
this diagram, historical service event data and the performance historical data are stored in a database at 100. After 
identifying a case study, the historical service event data and performance historical data are separated out according 
to the case study at 1 02. If the historical service event data and performance historical data are not in the case study, 
then this data are used to build a model as described in Fig. 2 at 104. Project service incidents along with statistical 

45 confidence bounds that should take place in the case study are determined at 1 06. An example of the project service 
incidents along with statistical confidence bounds are shown at 108. The service incidents are compared to data that 
are used in the case study at 11 0. If the actual data do not match the projection, then the model needs to be reexamined 
as noted at 112. On the other hand, if the data do match the projection within the statistical confidence bounds, then 
the model is validated as noted at 114. 

so [0036] The foregoing flow charts of this disclosure show the architecture, functionality, and operation of a possible 
implementation of the system for predicting the timing of future service events of a product. In this regard, each block 
represents a module, segment, or portion of code, which comprises one or more executable instructions for imple- 
menting the specified logical function(s). It should also be noted that in some alternative implementations, the functions 
noted in the blocks may occur out of the order noted in the figures, or for example, may in fact be executed substantially 

55 concurrently or in the reverse order, depending upon the functionality involved. 

[0037] The above-described system and method for predicting the timing of future service events of a product com- 
prise an ordered listing of executable instructions for implementing logical functions. The ordered listing can be em- 
bodied in any computer-readable medium for use by or in connection with a computer-based system that can retrieve 
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the instructions and execute them. In the context of this application, the computer-readable medium can be any means 
that can contain, store, communicate, propagate, transmit or transport the instructions. The computer readable medium 
can be, for example but not limited to, an electronic, magnetic, optical, electromagnetic, infrared system, apparatus, 
or device. An illustrative, but non -exhaustive list of computer-readable mediums can include an electrical connection 

5 (electronic) having one or more wires, a portable computer diskette (magnetic), a random access memory (RAM) 
(magnetic), a read-only memory (ROM) (magnetic), an erasable programmable read-only memory (EPROM or Flash 
memory) (magnetic), an optical fiber (optical), and a portable compact disc read-only memory (CDROM) (optical). It 
is even possible to use paper or another suitable medium upon which the instructions are printed. For instance, the 
instructions can be electronically captured via optical scanning of the paper or other medium, then compiled, interpreted 

10 or otherwise processed in a suitable manner if necessary, and then stored in a computer memory. 

[0038] For completeness, various aspects of the invention are set out in the following numbered clauses: 

I . A system for predicting the timing of a future service event of a product formed from a plurality of compartments, 
comprising: 

15 

a database that contains a plurality of service information and a plurality of performance information for the 
. product; 

a statistical analyzer that analyzes the plurality of service information to determine a plurality of compartment 
failure information; 

20 a performance deterioration rate analyzer that analyzes performance deterioration rate of the product from 

the plurality of service information and performance information; and 

a simulator for simulating a distribution of future service events of the product according to the plurality of 
compartment failure information and the performance deterioration rate analysis. 

25 2. The system according to clause 1, wherein the database comprises a service database and a performance 

historical database. 

3. The system according to clause 1 , wherein the plurality of service information comprises compartment defini- 
tions, repair history and service factors. 

30 

4. The system according to clause 1, wherein the plurality of performance information comprises performance 
characteristic values, initial data levels after servicing, current data levels, dates at which the product is serviced, 
and variables that affect the servicing of a subset of the plurality of compartments. 

35 5. The system according to clause 1, further comprising a preprocessor for processing the plurality of service 

information into a predetermined format. 

6. The system according to clause 1 , wherein the preprocessor generates a plurality of data files according to the 
plurality of service information. 

40 

7. The system according to clause 1 , wherein the plurality of compartment failure information comprises compart- 
ment failure variables and compartment time-to-failure coefficients. 

8. The system according to clause 7, wherein the statistical analyzer uses the plurality of compartment failure 
45 information to determine which compartment failure variables influence the timing of future service events and 

estimate time-to-failure distributions for the plurality of compartments. 

9. The system according to clause 8, wherein the statistical analyzer uses the estimated time-to-failure distributions 
to determine a Weibull distribution for a subset of the plurality of compartments defined for the product, 

10. The system according to clause 1, wherein the statistical analyzer comprises a service analysis script that 
executes a plurality of statistical procedures. 

II. The system according to clause 10, wherein the plurality of statistical procedures comprise a multivariate 
55 regression and/or a correlation analysis. 

12. The system according to clause 10, wherein the service analysis script generates a plurality of statistical di- 
agnostic information. 
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13. The system according to clause 12, wherein the plurality of statistical diagnostic information comprises good- 
ness-of-fit metrics and collinearity diagnostics. 

14. The system according to ciause 10, wherein the service analysis script generates a plurality of residual plots. 

1 5. The system according to clause 1 , wherein the statistical analyzer comprises a validation script. 

1 6. The system according to clause 1 5, wherein the validation script is applied to a plurality of case studies set up 
for the product. 

1 7. The system according to clause 1 , wherein the performance deterioration rate analyzer comprises a statistical 
analysis script that relates a subset of compartments of the product according to time. 

1 8. The system according to clause 1 7, wherein the statistical analysis script generates an estimated deterioration 
rate curve for the subset of compartments of the product. 

19. The system according to clause 18, wherein the performance deterioration rate analyzer further comprises a 
transformer that transforms each estimated deterioration rate curve for a compartment to a performance life dis- 
tribution. 

20. The system according to clause 19, wherein the simulator uses the performance life distributions to determine 
a Weibull distribution for a subset of the plurality of compartments defined for the product. 

21 . The system according to clause 1 , wherein the simulator forecasts a service plan for the future service events 
that comprises the time for scheduling the service events. 

22. A system for predicting the timing of a future service event of a product formed from a plurality of compartments, 
comprising: 

means for containing a plurality of service information and a plurality of performance information forthe product; 
means for analyzing the plurality of service information to determine a plurality of compartment failure infor- 
mation; 

means for performing a deterioration rate analysis that determines performance deterioration rate of the prod- 
uct from the plurality of service information and performance information; and 

means for simulating a distribution of future service events of the product according to the plurality of com- 
partment failure information and the performance deterioration rate analysis. 

23. The system according to clause 22, wherein the plurality of service information comprises compartment defi- 
nitions, repair history and service factors. 

24. The system according to clause 22, wherein the plurality of performance information comprises performance 
characteristic values, initial data levels after servicing, current data levels, dates at which the product is serviced, 
and variables that affect the servicing of a subset of the plurality of compartments of the product; 

25. The system according to clause 22, further comprising nieans for preprocessing the plurality of service infor- 
mation into a predetermined format. 

26. The system according to clause 25, wherein the preprocessing means generates a. plurality of data files ac- 
cording to the plurality of service information. 

27. The system according to clause 22, wherein the plurality of compartment failure information comprises com- 
partment failure variables and compartment time-to-failure coefficients. 

28. The system according to clause 27, wherein the analyzing means uses the plurality of compartment failure 
information to determine which compartment failure variables influence the timing of future service events and 
estimate time-to-failure distributions for the compartments. 

29. The system according to clause 28, wherein the analyzing means uses the estimated time-to-failure distribu- 
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tions to determine a Weibull distribution for a subset of the plurality of compartments defined for the product. 

» 

30. The system according to clause 22, wherein the analyzing means comprises a service analysis script that 
executes a plurality of statistical procedures. 

31. The system according to clause 30, wherein the plurality of statistical procedures comprise a multivariate 
regression and/or a correlation analysis. 

32. The system according to clause 30, wherein the service analysis script generates a plurality of statistical di- 
agnostic information. 

33. The system according to clause 32, wherein the plurality of statistical diagnostic information comprises good- 
ness-of-fit metrics and collinearity diagnostics. 

34. The system according to clause 30, wherein the service analysis script generates a plurality of residual plots. 

35. The system according to clause 22, wherein the analyzing means comprises a validation script. 

36. The system according to clause 35, wherein the validation script is applied to a plurality of case studies set up 
for the product. 

37. The system according to clause 22, wherein the performing means comprises a statistical analysis script that 
relates a subset of the plurality of compartments of the product according to time. 

38. The system according to clause 37, wherein the statistical analysis script generates an estimated deterioration 
rate curve for a subset of the plurality of compartments of the product. 

39. The system according to clause 38, wherein the performing means further comprises means for transforming 
each estimated deterioration rate curve for a compartment to a performance life distribution. 

40. The system according to clause 39, wherein the simulator uses the performance life distribution to determine 
a Weibull distribution for a subset of the plurality of compartments defined for the product. 

41 . The system according to clause 22, wherein the simulator forecasts a service plan for the future service events 
that comprises the time for scheduling the service events. 

42. A method for predicting the timing of a future service event of a product formed from a plurality of compartments, 
comprising; 

storing a plurality of service information and a plurality of performance information for the product; 
analyzing the plurality of service information to determine a plurality of compartment failure information; 
performing a deterioration rate analysis of the product from the plurality of service information and performance 
information; and 

simulating a distribution of future service events of the product according to the plurality of compartment failure 
information and the deterioration rate analysis. 

43. The method according to clause 42, wherein the plurality of service information comprises compartment def- 
initions, repair history and service factors. 

44. The method according to clause 42, wherein the plurality of performance information comprises performance 
characteristic values, initial data levels after servicing, current data levels, dates at which the product is serviced, 
and variables that affect the servicing of a subset of the plurality of compartments of the product, 

45. The method according to clause 42, further comprising preprocessing the plurality of service information into 
a predetermined format. 

46. The method according to clause 45, wherein the preprocessing generates a plurality of data files according to 
the plurality of service information. 
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47. The method according to clause 42, wherein the plurality of compartment failure information comprises com- 
partment failure variables and compartment time-to-failure coefficients. 

48. The method according to clause 47, wherein the analyzing uses the plurality of compartment failure information 
to determine which compartment failure variables influence the timing of future service events and estimate time- 
to-failure distributions for the plurality of compartments. 

49. The method according to clause 48, wherein the analyzing uses the estimated time-to-failure distributions to 
determine a Weibull distribution for a subset of the plurality of compartments. 

50. The method according to clause 42, wherein the analyzing comprises using a service analysis script that 
executes a plurality of statistical procedures. 

51. The method according to clause 50, wherein the plurality of statistical procedures comprise a multivariate 
regression and/or a correlation analysis. 

52. The method according to clause 51, wherein the service analysis script generates a plurality of statistical 
diagnostic information. 

53. The method according to clause 52, wherein the plurality of statistical diagnostic information comprises good- 
ness-of-fit metrics and col linearity diagnostics. 

54. The method according to clause 50, wherein the service analysis script generates a plurality of residual plots. 

55. The method according to clause 42, wherein the analyzing comprises using a validation script. 

56. The method according to clause 55, wherein the validation script is applied to a plurality of case studies set 
up for the product. 

57. The method according to clause 42, wherein the performing comprises using a statistical analysis script that 
relates a subset of the plurality of compartments of the product according to time. 

58. The method according to clause 57, wherein the statistical analysis script generates an estimated deterioration 
rate curve for a subset of the plurality of compartments of the product. 

59. The method according to clause 58, wherein the performing a deterioration rate analysis further comprises 
transforming each estimated deterioration rate curve for a compartment to a performance life distribution. 

60. The method according to clause 59, wherein the simulating uses the performance life distributions to determine 
a Weibull distribution for a subset of the plurality of compartments. 

61 . The method according to clause 42, wherein the simulating forecasts a service plan for the future service 
events that comprises the time for scheduling the service events. 

62. A computer-readable medium storing computer instructions for instructing a computer system to predict the 
timing of a future service event of a product formed from a plurality of compartments, the computer instructions 
comprising: 

storing a plurality of service information and a plurality of performance information for the product; 
analyzing the plurality of service information to determine a plurality of compartment failure information; 
performing a deterioration rate analysis of the product from the plurality of service information and performance 
information; and 

simulating a distribution of future service events of the product according to the plurality of compartment failure 
information and the deterioration rate analysis. 

63. The computer-readable medium according to clause 62, wherein the plurality of service information comprises 
compartment definitions, repair history and service factors. 
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64. The computer-readable medium according to clause 62, wherein the plurality of performance information com- 
prises performance characteristic values, initial data levels after servicing, current data levels, dates at which the 
product is serviced , and variables that affect the servicing of a subset of the plu rality of compartments of the product. 

65. The computer-readable medium according to clause 62, further comprising preprocessing instructions that 
preprocess the plurality of service information into a predetermined format. 

66. The computer-readable medium according to clause 65, wherein the preprocessing instructions generates a 
plurality of data files according to the plurality of service information. 

67. The computer-readable medium according to clause 62, wherein the plurality of compartment failure information 
comprises compartment failure variables and compartment time-to-failure coefficients. 

68. The computer-readable medium according to clause 67, wherein the analyzing instructions uses the plurality 
of compartment failure information to determine which compartment failure variables influence the timing of future 
service events and estimate time-to-failure distributions for the plurality of compartments. 

69. The computer-readable medium according to clause 68, wherein the analyzing instructions use the estimated 
time-to-failure distributions to determine a Weibull distribution for a subset of the plurality of compartments. 

70. The computer-readable medium according to clause 62, wherein the analyzing instructions comprises instruc- 
tions for using a service analysis script that executes a plurality of statistical procedures. 

71 . The computer-readable medium according to clause 70, wherein the plurality of statistical procedures comprise 
a multivariate regression and/or a correlation analysis. 

72. The computer-readable medium according to clause 71 , wherein the service analysis script generates a plurality 
of statistical diagnostic information. 

73. The computer-readable medium according to clause 72, wherein the plurality of statistical diagnostic informa- 
tion comprises goodness-of-fit metrics and collinearity diagnostics. 

74. The computer-readable medium according to clause 70, wherein the service analysis script generates a plurality 
of residual plots. 

75. The computer-readable medium according to clause 62, wherein the analyzing instructions comprise using a 
validation script. 

76. The computer-readable medium according to clause 75, wherein the validation script is applied to a plurality 
of case studies set up for the product. 

77. The computer-readable medium according to clause 62, wherein the performing instructions comprise using 
a statistical analysis script that relates a subset of the plurality of compartments of the product according to time. 

78. The computer-readable medium according to clause 77, wherein the statistical analysis script generates an 
estimated deterioration rate curve for a subset of the plurality of compartments of the product. 

79. The computer-readable medium according to clause 78, wherein the performing instructions further comprise 
transforming instructions that transform each estimated deterioration rate curve to a performance life distribution. 

80. The computer-readable medium according to clause 79, wherein the simulating instructions use the perform- 
ance life distribution to determine a Weibull distribution for a subset of the plurality of compartments. 

81 . The computer-readable medium according to clause 62, wherein the simulating instructions forecasts a service 
plan for the future service events that comprises the time for scheduling the service events. 
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Claims 

I . A system for predicting the timing of a future service event of a product formed from a plurality of compartments, 
comprising: 

5 

a database that contains a plurality of service information and a plurality of performance information for the 
product; 

a statistical analyzer that analyzes the plurality of service information to determine a plurality of compartment 
failure information; 

10 a performance deterioration rate analyzer that analyzes performance deterioration rate of the product from 

the plurality of service information and performance information; and 

a simulator for simulating a distribution of future service events of the product according to the plurality of 
compartment failure information and the performance deterioration rate analysis. 

15 2. A method for predicting the timing of a future service event of a product formed from a plurality of compartments, 
comprising; 

storing a plurality of service information and a plurality of performance information for the product; 
analyzing the plurality of service information to determine a plurality of compartment failure information; 
20 performing a deterioration rate analysis of the product from the plurality of service information and performance 

information; and 

simulating a distribution of future service events of the product according to the plurality of compartment failure 
information and the deterioration rate analysis. 

25 3. The method according to claim 2, wherein the plurality of service information comprises compartment definitions, 
repair history and service factors. 

4. The method according to claim 2, wherein the plurality of performance information comprises performance char- 
acteristic values, initial data levels after servicing, current data levels, dates at which the product is serviced, and 

30 variables that affect the servicing of a subset of the plurality of compartments of the product. 

5. The method according to claim 2, further comprising preprocessing the plurality of service information into a pre- 
determined format. 

35 6. The method according to claim 2, wherein the plurality of compartment failure information comprises compartment 
failure variables and compartment time-to-failure coefficients. 

7. The method according to claim 2, wherein the analyzing comprises using a service analysis script that executes 
a plurality of statistical procedures. 

40 

8. The method according to claim 2, wherein the analyzing comprises using a validation script. 

9. The method according to claim 2, wherein the performing comprises using a statistical analysis script that relates 
a subset of the plurality of compartments of the product according to time. 

45 

10. The method according to claim 2, wherein the simulating forecasts a service plan for the future service events that 
comprises the time for scheduling the service events. 

I I . A computer-readable medium storing computer instructions for instructing a computer system to predict the timing 
50 of a future service event of a product formed from a plurality of compartments, the computer instructions comprising: 

storing a plurality of service information and a plurality of performance information for the product; 
analyzing the plurality of service information to determine a plurality of compartment failure information; 
performing a deterioration rate analysis of the product from the plurality of service information and performance 
55 information; and 

simulating a distribution of future service events of the product according to the plurality of compartment failure 
information and the deterioration rate analysis. 
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FIG. 2 
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FIG. 3 
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FIG. 4 
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FIG. 5 
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